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We discuss the possibility of a Kondo like effect associated with H in metals resulting from the 
strong dependence of the H Is orbital radius on the occupation number. We demonstrate that such 
a strong breathing property of the orbital radius, which translates directly into a strong occupation 
dependent hopping, results in the formation of local singlet-like bound states involving one electron 
on H and one on the surrounding metal orbitals. We also show that already at a mean field 
level an occupation dependent hopping integral leads to a substantial potential energy correction on 
hydrogen, and that the failure of band structure methods to incorporate this correction is responsible 
for the incorrect prediction of a metallic ground state for the YH3 switchable mirror compounds. 



I. INTRODUCTION 

The recent discovery j^] of the so called switchable mir- 
ror compounds based on YH^-x has renevifed the interest 
in the electronic structure of transition metal and rare 
earth hydrides. These compounds undergo a metal to 
insulator transition as x changes from one to zero with 
the accompanying change from high optical reflectivity 
to optical transparency for photon energies below about 
2 eV. LDA band structure calculations fail to reproduce 
the semiconducting gap for a; = 3 in the widely accepted 
H0D3 structure or in the simpler cubic LaHs structure. 
Although a gap can be obtained for more complicated 
distorted structures it is still much too small; and 
although such distortions are not excluded for YH3, they 
do not appear to occur for LaHa, which has similar prop- 
erties. 

The failure to produce large enough band gaps in semi- 
conductors by LDA is a well-known shortcoming which 
here however seems to take quite dramatic forms: The 
valence and conduction bands in LDA overlap by about 
1 eV, so a total relative shift of about 3 eV is required 
to match the experimental value. Such dramatic discrep- 
ancies are reminiscent of the strongly correlated systems 
like the transition metal and rare earth oxides and point 
perhaps to the importance of correlation effects. Since 
however the H Is orbitals are rather extended as com- 
pared to the id's, of the transition metals, especially for 
the negative ion, they are expected to form rather broad 
bands, and the on-site Coulomb interactions are strongly 
screened. Therefore the origin of the correlation effects 
may be quite different. In this paper we address this 
problem and come to the new suggestion that the cor- 
relation effects are a consequence of the large change in 
the H Is orbital radius upon orbital occupation. This 
"breathing" property of the hydrogen ion is shown to in- 
troduce a new term in the mean field treatment of the 
electronic structure of hydrides. This term results in an 
opening up of the band gap in a quite natural way, with 
the retention of large band widths and nearly one particle 



behaviour of the excited states. Using a model Hamilto- 
nian to demonstrate this behaviour we also show that for 
a range of parameter values the system behaves like a 
Kondo lattice insulator similar to that suggested by Ng 
et al 101. 



II. THE BREATHING HYDROGEN ATOM 

As is well known and referred to in most general chem- 
istry text books, the so called effective radius of hydrogen 
is extremely strongly dependent on the charge state. The 
crystal radius of neutral H is 0.26 A, whereas that of the 
negative ion H~ is 1.54 A The values of the aver- 
age Is orbital radius \/<r^ >, as obtained from free ion 
Hartree-Fock calculations, are 0.8 A for Hand 1.72 A for 
. This very large change is not unexpected since in 
H, with its low nuclear charge of one, the screening of the 
nuclear Coulomb potential by a second s electron is very 
important and has a dramatic effect on the orbital radius. 
This large effect causes the effective hopping integrals or 
hybridizations with surrounding ions to be strongly dif- 
ferent for the fluctuations involving H to H+ as compared 
to those involving H to H^, as pictured in Fig. 0. 
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FIG. 1. Schematic representation of two electron and single 
electron wave function with (left column) and without (right 
column) taking into account the expansion of the Hydrogen 
wave function. 
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If these instantaneous changes in the hopping integrals 
are larger than or comparable to other energy scales, like 
the orbital energy splittings, they must be treated explic- 
itly. They can not be treated in a mean field like way (as 
is done in band theory), taking an average orbital radius 
corresponding to the average occupation as determined 
from self-consistent calculations. 

To model the effect of the "breathing" hydrogen we in- 
troduce an occupation dependent hopping integral (fol- 
lowing Hirsch between H and its nearest neighbours, 
in addition to the usual on-site Coulomb interaction of 
the Hubbard or Anderson impurity |^ models. We 
consider the Hamiltonian of the form: 

= -A ^ hlK + Uh\h] hlhi 

a 

- ^ [ ( VillKhlhs + V^llKhshl ) + H.c. ], (1) 

a 

which describes the hybridization of a single hydrogen 
atom with a single ligand orbital. Here hi (Zj) creates an 
electron on the hydrogen (ligand), there is a charge trans- 
fer energy A between hydrogen and ligand (we assume 
A > 0) and a Coulomb repulsion U between electrons 
on hydrogen. Vi and V2 are the occupation dependent 
hydrogen-ligand hopping integrals. We discuss the first 
ionization energy E- for photoemission on a single cell, 
consisting of a hydrogen atom and a ligand orbital, filled 
with two electrons. The minimum ionization energy 
together with the minimal electron affinity E+ , determine 
the excitation gap Egap — -E- — -E+ . An underestimation 
of the magnitude of E- thus may lead to a too small gap 
energy, as seems to be the case in LDA calculation for 
YH3. 

For simplicity we take V2 = 0, so that the ground state 
of a single electron just corresponds to the electron be- 
ing trapped in the collapsed hydrogen orbital, and has 
energy —A. A straightforward calculation then gives the 
corresponding ionization energy 
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Taking for simplicity [/ = A we obtain E = -^/2Vi. The 
ionization energy thus is predominantly due to the loss 
of kinetic energy, because the single electron in the final 
state cannot escape from the collapsed hydrogen orbital, 
so that the large gain in kinetic energy, which was possi- 
ble for two electrons, is no longer possible. 
On the other hand, constructing a single-particle Hamil- 
tonian with an averaged hopping integral Vmf ~ 
{h1hcr)Vi and an "effective" on-site energy Amf (as it is 
done in an LDA calculation), the excitation energy would 
be simply the energy of the occupied mean-field orbital. 
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If the occupation of hydrogen Is is significantly smaller 
than one per spin state, this way of calculation will miss 
a large part of the kinetic energy contribution to the ex- 
citation energy, unless the "effective" on-site energy is 
corrected to take this effect into account. 
In a mean-field treatment of this Hamiltonian, it will be- 
come apparent that the occupation dependent hopping 
gives rise to very peculiar physics. Breaking down the 
conditional hopping terms into quadratic terms we ob- 
tain: l^hahlhg l\ha{h\hg) + {l1h„)h\hg. The first 
of these terms corresponds to weighting the "large" hy- 
bridization integral by the occupation of the hydrogen 
orbital, which is what one might have expected; the sec- 
ond term, however, is a correction to the on-site energy 
of hydrogen by a part of the kinetic energy. All in all we 
obtain: 

Hmf = [ -^MFhlK + ( VmfIIK + H.c. ) ], 



Vmf = VinH + V2 (1 - nn) 
^Amf = -A + U uh + a{T) 
V1-V2 



a 



if 



(2) 



where ~ {h\h^) and (T) is the energy of hybridiza- 
tion between hydrogen and ligand. We thus find the sur- 
prising result that in this approximation the expectation 
value of the kinetic energy (T) enters as an additional 
"potential" on the hydrogen sites, and is in fact even 
enhanced by the factor a. In the limit Vi 3> V2 we find 
a — !■ n'^ > 1, so that the correction to the on-site poten- 
tial of hydrogen becomes {Vi/Vmf){T) , i.e., the kinetic 
energy for mixing with the hydrogen site, but calculated 
with the hopping integral for the "expanded" atom. This 
is clearly a huge energy, but it has a very clear physical 
interpretation: for two electrons in the cell, the hydro- 
gen atom will oscillate between and H~ , so as to 
take maximum advantage of the expansion of the wave 
function, and the hybridization energy will be large. Re- 
moving one electron, the remaining electron will essen- 
tially be trapped in the "collapsed" hydrogen orbital, and 
there is practially no more hybridization energy. In the 
mean-field wave function both electrons are in the lower 
molecular orbital, which (due to its strongly negative ef- 
fective on-site energy resulting from Eq. ^ has predom- 
inant hydrogen character. The ionization energy, which 
by Koopmans theorem should be given by the mean-field 
eigenvalue, then contains almost the entire kinetic energy 
of the two-electron state, because this kinetic energy has 
been put into the on-site potential of the hydrogen atom. 
To make this more quantitative, we have performed ex- 
act diagonalization calculations for a ID chain of a model 
with occupation dependent hybridization between hy- 
drogen and ligand. A schematical representation of the 
model is given in Fig. ^ 
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FIG. 2. Schematic representation of the ID model used in 
the exact diagonalization. 



Computer memory and CPU time limitations prohibit 
to diagonalize chains of more than six unit cells of this 
model, at least for the most interesting densities near 
two electrons per unit cell. For the given parameters, the 
ground state of this system corresponds closely to a state 
with two electrons in a local singlet state in each unit cell, 
with one electron primarily on H and the other primarily 
on the neighbouring ligands. Such a state is reminis- 
cent of a Kondo-lattice insulator ground state. Figure ^ 
shows the fc-resolved electron addition and removal spec- 
tra for the six-site chain at half-filling. To get a feeling for 
the dispersions of the different bands, we have combined 
spectra calculated with both periodic and antiperiodic 
boundary conditions. While there is no rigorous proof 
for this, inspection of Fig. ^ shows that in this way one 
obtains remarkably smooth "band structures" . 



Unlike spectra for, e.g., the Hubbard model at half-filling, 
the calculated spectra are surprisingly "coherent", with 
almost all of the spectral weight being concentrated in 
just three well defined "bands" . In inverse photoemis- 
sion there is a band of predominant metal character, with 
little dispersion. In photoemission, there is a strongly 
dispersive band of mixed hydrogen-metal character, and 
a dispersionless low intensity band of practically pure 
hydrogen character. More detailed analysis shows, that 
the dispersionless band corresponds to final states 
(i.e., it is a kind of "lower Hubbard band") whereas the 
dipersive band corresponds to iJ°-like final states. Next, 
Fig. ^ shows the spectral function for different values 
of U. For comparison, the bands obtained by a mean- 
field solution of the model are also shown. Thereby the 
calculation has been done both with and without the 
kinetic energy correction Vt — ct{T) to the hydrogen on- 
site potential. Quite obviously the calculation with Vt 
reproduces the exact band structure very well, whereas 
the bands without Vt while giving roughly the correct 
dispersion substantially underestimate the gap size. As 
explained above, we believe that LDA misses the kinetic 
energy correction Vt, so that the LDA band structure 
rather corresponds to the bands without Vt- In a phe- 
nomenological way, this suggests a kind of "scissors op- 
erator" approach to obtain the "correct" band structure 
of YH3 from the LDA result. 
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FIG. 3. Single particle spectral function of the ID model 
with six unit cells. The part to the right (left) of the thin ver- 
tical line correspond to electron addition (removal) from the 
half-filled ground state (i.e., two electrons per unit cell). The 
full line corresponds to electron removal or addition on the 
"metal" sites, the dashed line to hydrogen. Parameter values 
are Vx = 2 eV, Vi = 0.2 eV, t = 0.5 eV, A = 1 eV (see Fig. 

The full (dashed) dispersion curve gives the mean-field 
bands calculated with (without) the correction Vt = 



III. APPLICATION TO YH3 



We now use the above ideas for the case of YH3 and 
attempt to obtain reasonable parameters and subsequent 
estimates of the band gap. The three 5c? electrons of Y 
will in the above scenario all be bound by the three H 
atoms per Y; this would again result in an insulating 
ground state. First, we obtain good estimates for the 
average hopping integrals and on-site energies, using a 
tight binding fit to an LDA band structure calculated 
with the LMTO method The upper panel of Fig. | 
shows the band structure calculated for YH3 in the LaHa 
structure, which is practically identical to previous pub- 
lished results In the lower panel we show the 
result of a calculation for Y metal with the lattice con- 
stants of YH3, in order to establish the contribution of 
the H Is orbitals. 
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FIG. 4. Band structures YH3 (top) and metallic yttrium 
with the lattice constants of YH3 (bottom). The Fermi en- 
ergy is marked by the dashed line. 



To begin with, the lowest two bands which have pre- 
dominant 5sp character, are nearly identical in both 
solids. Obviously, these free-electron-like states are not 
significantly affected by the insertion of hydrogen into 
the interstitial sites. In yttrium metal, the next group of 
states are the Y Ad bands, which form a dense continuum 
with a width of ~ 6 eV. The Fermi energy cuts into the 
lower part of this complex of d bands. In YII3 the sit- 
uation is very different: the lowest Y 4d band (which is 
still quite reminiscent of the third-lowest band in metal- 
lic yttrium) is split off from the remaining d bands by 

4 eV, and in the resulting gap bands of predominantly 
hydrogen character are inserted. Although these hydro- 
gen bands do have an appreciable width, they barely 
overlap with the lowermost of the following d bands. In 
fact, the band structure of YH3 already shows a very 
clear "gap" between the top of the hydrogen like valence 
band and the Y Ad-like conduction band throughout the 
entire Brillouin zone - YII3 thus is already "almost" a 
semiconductor. The Fermi energy cuts into the top of 
the hydrogen-like valence band and the bottom of the d- 
like conduction band, so that LDA predicts YH3 to be a 
semimetal. The shift of the lowermost d to considerably 
higher binding energy upon insertion of hydrogen, which 
is predicted by LDA, is in qualitative agreement with 
the photoemission data of Fujimori and Schlapbach [p^ : 
For metallic yttrium, these authors found a high inten- 
sity structure at binding energies < 2 eV, which probably 
corresponds to the occupied part of the Y 4d bands. For 
YII3 a similar structure occurs at a binding energy of 
6 eV, indicating the shift of the d band away from Ep. 
It should be noted, however, that the experimental shift 
is larger by « 2 eV than that predicted by LDA. The 
picture thus is quite reminiscent of the well-known band- 
gap problem in semiconductors, where LDA fails to give 
correct values for the semiconducting gaps. 
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FIG. 5. LDA band structure for YH3 (full line) and tight 
binding fit (dashed line). Energies are in eV. 

To extract additional information, we performed a 
tight-binding fit to the LDA band structure. It turned 
out that by using hybridization integrals only between 
nearest neighbours a surprisingly good fit of the first 
few valence and conduction bands could be obtained, as 
shown in Fig. ^. The tight binding fit was obtained 
using the following Hamiltonian in a mean field way, as 
discussed below: 

h.v^a j,(T 

+ Uj2h]^^hj,^h]fy^^. (3) 
j 

Here e,y{k) denotes the Y Ad bands, which we obtained 

from the tight-binding fit. 

In a mean-field treatment (^ would turn into 

H = Y1 e^(fc)4,^,<,dfe,z.,<T + ^3h],ah3,<y 

V^(S),/U/»..-(l-".)) + ^^-c.], (4) 

with rij = {hj -^hj^^) = {h}- ^hj^i) and Aj = ij -\- Urij. 
We now introduce the parameter A, which we assume in- 
dependent of j, as A = v/-V,^ J%^i^l\ o j i-e-, the ratio of 
hopping integrals for the collapsed and expanded hyb- 
drogen atom. We then estimate the change as 

(5) 

where V^"f j is the hybridization integral extracted from 
the tight-binding fit. Since rij can be obtained from the 
tight-binding calculation as well, we can thus, for given 
A, obtain an estimate of vj^^K . Next, we estimate the 
"bare" on-site energies ej of the hydrogen atoms from 
those of the tight-binding fit, Aj^, as follows: 
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(6) 



then is shown in Fig. 



This introduces another unknown parameter: the onsite 
Coulomb repulsion U between two electrons on the hy- 
drogen site. While an ab-initio calculation of U and A 
would be highly desirable, this is outside the range of 
techniques available to us. We therefore will treat these 
quantities as implicit parameters, and consider the vari- 
ation of possible results when U and A are varied within 
"reasonable bounds" . 

Using the parameters estimated in this way, we now pro- 
ceed to an impurity-like calculation to determine the sta- 
bilization energy of Kondo-like local singlets, formed on 
a single hydrogen atom in the lattice of Y 4c? orbitals. In 

(2) 

the first step, we calculate the ground state energy Eg ' 
of a two-electron bound state from the ansatz 

= [a/,t/,t + -Lg/?,( dl ,^^^hl + h]dl^^^ ) ]\vac). 

Here summation over k and v refers to the Y Ad bands. 
The dispersion of these bands and the hybridization ma- 
trix elements between the band states and the hydrogen 
atom are calculated using the parameters from the tight 
binding fit, whereby the hydrogen-yttrium hybridizations 
depend on A and the hydrogen on-site energy on U. 
Then, we want to know the stability of this state against 
decay into a state with a single electron remaining in 
the hydrogen atom, and the second electron being in a 
free yttrium c?-like state. The energy iJg^-' of the single 
electron in the hydrogen is calculated from the ansatz 

1*^)^ [a'h\+J2pldl^^^^]\vac), 

and for the energy of the d-like electron we simply choose 
the lower bound of the d-band complex, E^dge ■ We then 
form the difference AS = £;^^^ - {e'^'^ -t- E^dge) (see Fig. 




FIG. 6. Schematic representation of the stabilization en- 
ergy for the local two-electron bound state. 

which obviously determines the stability of the two- 
electron state against decay. This energy will be a func- 
tion of the unknown parameters U and A. The result 



One can see that for "reasonable" values of U and 
not very extreme values of A the two-electron bound 
state attains a stabilization energy of several electron- 
volts. Drawing an analogy with the situation in cuprate 
superconductors, where the Zhang-Rice singlet has a sta- 
bilization energy of approximately 1 eV, it seems quite 
reasonable to adopt the picture of local bound states. 
Then, for YH3 one may expect that these bound states 
form a split-off band, with the Fermi energy lying in the 
gap between these states and the bottom of the Ad band; 
the physics is similar to our exact result on the ID clus- 
ter. 
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FIG. 7. Ionization energies AE calculated with the impu- 
rity model for the metal plane hydrogen (left panel) and the 
tetragonal hydrogen (right panel). Note that the other tight 
binding parameters depend on the values for U and A through 
Equations (^ and 



IV. DOPING DEPENDENCE OF THE 
ELECTRONIC STRUCTURE 

To simulate the physics of what happens as one re- 
moves hydrogen from YH3, we resort to an exact cal- 
culation of a small cluster as done above, but we now 
remove one hydrogen site and one electron. In Fig. ^ we 
show the spectral function for electron removal and elec- 
tron addition before and after removal of the hydrogen. 
In the top panel, one can clearly recognize the large gap 
between the hydrogen like valence band and the metal 
like conduction band. Upon hydrogen removal, as shown 
in the bottom panel, the Fermi energy jumps into the 
metal band which implies that hydrogen behaves like an 
ion in that it binds two electrons. This is consistent 
with the above discussion and also suggested by Ng et al 

i- 

In our previous discussion we came to the conclusion 
that H actually binds two electrons: one localized on 
H, and the other on the nearest neighbour metal atoms. 
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With the removal of a H atom, which takes only one 
electron with it, another "lonely" electron is left behind, 
which must then be in the conduction band. Removal 
of H from the trihydride insulator should then transfer 
spectral weight for electron removal from the top of the 
valence band to the bottom of the conduction band, some 
2 eV higher in energy. This kind of behaviour upon dop- 
ing is very similar to that predicted for and observed 
iQ in the high Tc cuprates. Consistent with this are the 
observations by Peterman et al. ]l5[ |, who found that in 
hydrogen depleted trihydrides the Fermi energy falls into 
a "band" with very weak spectral weight, which grows 
upon further depletion. 

V. CONCLUSIONS 

We have shown that hydrogen is an extreme example 
of an atom with a large orbital occupation dependence 
of the orbital radius, leading to large occupation depen- 
dent hopping integrals in hydrides. We argue that ex- 
plicit inclusion of such terms in the Hamiltonian results 
in a "scissors operator" like separation of the valence and 
conduction bands and, consequently, the opening of a 
substantial gap. We argue that the insulating character 
of YH3 can be understood in this way. Using reasonable 
parameters obtained from tight binding fits to the band 
structure, we find that the ground state of YH3 corre- 
sponds closely to that of a Kondo insulator, with each H 
binding two electrons in a singlet state. 
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FIG. 8. Single particle fc-integrated spectral function for 
a three unit cell cluster of the ID model with open bound- 
ary conditions. The spectra are calculated at "half-filling" 
(upper part) and with one charge neutral hydrogen atom re- 
moved from the central cell (lower part). The parts of the 
spectra to the right (left) of the vertical dashed-dotted fine 
correspond to electron addition (removal). Parameter values 
are A = 1 eV, 17 = 2 eV, t = 0.5 eV, Vi = 2 eV, V2 = 0.2 eV. 
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